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Analytical  Models  for 
Battlespace  Information  Operations 

(BAT-IO) 

Donald  P.  Gaver 
Patricia  A.  Jacobs 

Part  2 

Battlespace  Information  War: 

Attack  and  Defense  of  a  Region  with  Adaptive 
Opponent  Behavior,  Including  Deterrence 

EXECUTIVE  SUMMARY 

The  coordination  of  information  acquisition  and  interpretation  to  direct  force 
application  is  increasingly  recognized  as  a  crucial  military  systems  design  and 
investment  issue.  This  paper  illustrates  tradeoffs  between  Blue/own  regional 
Attacker  sensor  and  shooter  capabilities:  it  studies  a  deep  strike  or  SCUD- 
hunting  scenario  in  a  low-resolution,  aggregated  manner  using  an  analytical 
state-space  approach  that  recognizes  gross  aggregated  regional  Defender  (Red), 
and  regional  Attacker  (Blue),  system  capabilities  and  limitations.  Emphasis  is 
accordingly  placed  on  explicitly  modeling  the  availability  and  utilization  of 
information  to  a  striking  Attacker,  as  it  becomes  available  from  a  realistically 
finite  sensor  and  C2  capability.  The  (imperfect)  information  on  opposition  units, 
the  Defenders,  that  are  candidates  for  prosecution  by  the  Attackers  is  passed  to 


the  finite,  hence  saturable  (here  missile-firing)  Attacking  force,  the  shooters,  that 
then  responds  by  prosecuting  those  units. 

The  models  specifically  recognize  that  regional  Defenders  will  not  be  detected 
immediately,  nor  recognized  perfectly,  nor  are  Defender  shots  (e.g.  SCUD 
launcher)  fired  perfectly,  or  immediately.  Furthermore,  attempts  to  effectively 
target  are  also  realistically  modeled  as  afflicted  by  imperfect  Attacker  battle  damage 
assessment  (BDA),  an  incapability  that,  if  pronounced,  will  tend  to  non-linearly 
saturate  shooters,  increase  their  response  times,  and  hence  reduce  targeting 
effectiveness  and  efficient  ammunition  expenditures.  Such  models  can  allow  for 
adaptation  by  both  attackers  and  defenders  to  recent  fortunes:  if  Defender 
presence  and  activity  is  effectively  countered  by  Attackers,  then  the  former  may 
tend  to  be  deterred  or  withdraw;  if  not,  the  Defenders  are  motivated  to  press 
their  apparent  advantage.  Sharp,  threshold-like,  responses  can  follow  from  the 
possibly  multi-stable  dynamics.  This  behavior  will  be  explored  in  this  report. 

The  present  models  are  mainly  deterministic  or  pseudo-stochastic  in  that  they 
represent  the  non-linear  ejfect  of  stochastic  saturation  approximately,  but 
adequately.  However,  they  can  straightforwardly  be  "made  stochastic", 
especially  Markovian,  and  so  realized  using  Monte  Carlo  simulations.  Computer 
programs  exist  to  provide  numerical  results;  some  are  given.  A  simple  one¬ 
dimensional  stochastic  (Markov  birth-death)  model  is  given  as  an  appendix  to 
this  paper.  This  model  can  be  shown  analytically  and  numerically  to  exhibit 
"stochastic  bi-stability"  properties  (two  attractors)  that  under  certain  circum¬ 
stances  (parameter  combinations)  lead  to  bimodal  steady-state  distributions  (if 
such  are  allowed  to  happen  by  the  dynamics,  and  are  of  interest).  Such  a 
tendency  will  occur  also  in  more  detailed,  but  less  analytically  tractable  models.  < 

There  are  many  problem  elements  that  have  been  initially  and  purposefully 
ignored.  They  are  addressed  in  later  work;  see  Gaver  and  Jacobs  (1999).  For 
instance,  the  effects  of  different  target  t5q3es,  false  targets,  and  decoys  must  be 
added  (some  "decoys"  are  in  effect  present,  in  the  form  of  killed  Defenders,  not 
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so  recognized,  that  are  mistakenly  re- targeted).  The  effect  of  different  principles 
for  Attacker  target  prioritization  under  uncertainty,  i.e.  dynamic  scheduling, 
requires  systematic  attention.  In  the  present  models  Attackers  are  invulnerable  to 
attack;  this  is  not  always  realistic,  and  can  be  changed  to  a  duel-like  scenario 
involving  suppression  of  enemy  (Defender)  air  defense  (SEAD);  a  first  paper  on 
this  topic  is  Gaver  and  Jacobs  (1998).  In  the  current  paper  Attackers  employ 
generic  missiles  only,  but  the  use  of  (vulnerable,  manned)  Attack  aircraft  can 
similarly  be  modeled,  as  can  combinations  of  Attack  aircraft.  Naval  gunfire,  and 
missiles,  recognizing  the  coordination  difficulties.  Employment  of  cued 
reconnaissance  aircraft,  possibly  UAVs,  can  likewise  be  represented 
quantitatively  as  state-space  components.  In  addition,  refinements  that  more 
faithfully  represent  spatial  and  perhaps  other  environmental  constraints  can  be 
incorporated,  as  can  details  of  communications  assets  and  message-handling 
protocols  in  use  by  both  Attackers  and  Defenders. 

The  present  paper  describes  some  of  the  possibilities  for  insights  inherent  in 
an  enhanced  state-space  approach.  As  pointed  out,  many  elaborations  are 
possible.  The  objective  is  to  recognize  only  that  detail  in  the  (preliminary)  models 
that  is  sufficient  to  hint  at  payoff  from  adding  suitable  assets  and  strategies  at 
appropriate  points  in  the  entire  system.  Finer  detail  and  resolution  is  left  to 
others  to  include,  and  possibly  profit  by.  More  elaborate  and  high-resolution 
models  within  such  tools  as  NSS  (METRON),  and  JWARS  eventually  can  focus 
with  greater  intensity  on  some  of  the  issues  raised  here. 

In  general  we  believe  that  this  report  is  in  accord  with  many  of  the  views  and 
suggestions  of  Ilachinski  (1996),  and  also  of  Dockery  and  Woodcock  (1993),  and 
others.  Those  two  publications  coritain  many  references,  some  to  previous  work 
on  related  topics.  Our  emphasis  on  explicit  representation  of  protagonists' 
information  state  is  of  interest  at  the  time  of  writing. 


Battlespace  Information  War: 

Attack  and  Defense  of  a  Region  with  Adaptive 
Opponent  Behavior,  Including  Deterrence 

(BAT-IW) 

Part  2 

D.P.  Gaver 
P.A.  Jacobs 
M.A.  Youngren 

1.  Introduction:  Generic  Adaptive  Staging 

Suppose  a  particular  territorial  region,  %  is  selected  by  an  attacking  force, 
denoted  as  Red,  as  one  possible  stage  from  which  to  assemble,  maneuver,  and 
possibly  launch  attacks.  An  example  is  the  designation  of  as  a  region  in  Iraq 
from  which  TELS/SCUDS  missile  systems  laimched  attacks  on  Israel  at  the  time 
of  the  Gulf  War;  another  such  region  was  the  staging  arena  for  attacks  on  Saudi 
Arabia.  But  the  region  could  equally  well  be  occupied  by  various  armed  units 
that  would  converge  on  a  particular  location  within  or  on  the  boundary  of  ^in 
case  it  were  attacked,  e.g.  by  amphibious  landing  operation. 

Let  the  force  opposing  Red  utilization  of  region  <^be  denoted  generically  as 
Blue.  In  the  present  scenario  Blue's  assets  consist  of  C4ISR  systems  ("sensors" 
plus  a  communication  system)  capable  of  detecting,  tracking  and  targeting  Red 
assets,  plus  weapons  for  attrition  (or  suppression)  of  those  assets  (the 


"shooters").  Importantly  and  realistically,  the  Blue  sensors  and  shooters  are 

I 

represented  as  of  limited  and  fallible  capability:  Red  forces  in  <^are  detected  and 
targeted  by  Blue  C4ISR  only  after  delay,  and  then  with  occasional  error,  and 
shooter  assets  are  sometimes  saturated  and  hit  and  kill  Red  assets  with  less  than 
perfect  success.  Shooter  assets  considered  here  are  a  generic  surface-surface 
missile  system,  e.g.  ATACMS;  with  some  additional  effort  a  force  of  manned 
aircraft  could  alternatively,  or  as  a  complement,  be  envisioned  and  modeled. 
Included  in  the  above  model  properties  is  the  explicit  recognition  that  follow-up 
of  a  Blue  attack  on  a  detected  Red,  i.e.  battle  damage  assessment  (BDA),  is  error- 
prone:  not  only  may  a  dead  target  be  mis-identified  as  alive  and  wastefully 
retargeted,  acting  effectively  as  a  decoy,  but  an  alive  target,  the  subject  of  an 
unsuccessful  attack,  can  be  ignored  and  hence  allowed  to  profitably  reposition  or 
re-attack. , 

1.1  Adaptive  Staging 

It  is  natural  to  expect  that  the  conflict  scenario  portrayed  occur  dynamically. 
For  example.  Red  forces  may  first  infiltrate  ^with  intention  of  establishing  a 
prescribed  strength  there,  perhaps  occasionally  carrying  out  individual  attacks. 

In  response.  Blue  forces  will  focus  increased  sensor  and  shooter  attention  on  the 

« 

region.  If  such  resistance  is  quick  and  effective  Red  activity  will  tend  to  adapt  by 
tactics  that  reduce  vulnerability,  but  also  effectiveness:  they  will  be  deterred  (at 
least  in  the  region  imder  consideration).  Correspondingly,  Blue  force  attention  to 
the  region  may  well  then  be  limited,  perhaps  re-directed.  Red  occupancy  of  ^ 
settles  at  an  annoyance  level.  On  the  other  hand,  if  Red  forces  once  surpass  the 
combined  sensor-shooter  forces  of  Blue  the  region  will  tend  to  settle  at  a 
considerably  higher  threat  level  -  one  considerably  more  difficult  for  Blue  to 
dislodge. 
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The  purpose  of  this  paper  is  to  connect  various  aspects  of  the  above 
information-influenced  scenario  described  by  simple  dynamic  models  that 
dramatically  exhibit  the  sensitivity  of  conflict  outcome  to  opponent's  assets  and 
operational  behavior.  The  models  proposed  for  first  study  are  aggregated  and 
deterministic,  but  still  provocative.  They  can  both  be  readily,  and  informatively, 
extended  to  far  more  detailed  higher  resolution  tools  that  give  more  refined  and 
hence  perhaps  useful  (or  at  least  easily  defensible)  guides  to  investment  decisions 
(e.g.  sensor  and  weapons  mix,  target  prioritization,  etc.).  The  approaches  of  the 
present  paper  provide  a  quick  preliminary  appraisal  of  a  hypothetical  situation; 
polish-up  can  be  accomplished  subsequently. 

2.  Model  I:  The  "SCUDWORM" Analogy 

About  twenty  years  ago,  D.  Ludwig  identified  and  explored  an  intriguing 
dynamic  situation  pertaining  to  the  interlinked  behavior  of  a  pest,  the  spruce 
budworm,  and  birds  that  feed  upon  populations  of  them.  See  Ludwig  et  al.  (1978) 
and  more  recently  Murray  (1989).  When  spruce  are  hospitable,  presumably  when 
buds  appear,  populations  of  budworms  appear  to  feed  upon  them  and  grow;  in 
turn,  the  population  of  worms  attracts  birds  that  consume  the  worms,  with 
attractiveness  increasing  with  worm  population  size.  If  both  budworm  and  bird 
population  sizes  are  ultimately  limited  the  dynamics  are  sometimes  consistent 
with  two  (non-zero)  local  stable  points  for  worm  (and  local  bird)  populations:  one 
low,  and  one  high.  The  simple  deterministic  model  implies  that  once  (initial) 
conditions  bring  worm  and  bird  populations  near  either  such  point  they  remain 
there  indefinitely,  even  though  all  descriptive  rate  parameters  are  otherwise 
unchanged.  Only  by  adding  more  to  the  problem,  e.g.  by  modifying  the  process 
to  become  stochastic  or  by  incorporating  time  dependence,  will  escape  from 


3 


quasi-steady  state  occur;  this  may  occur  as  a  result  of  smallish  change  of 
conditions,  and  hence  appear  instantaneous. 

We  exploit  the  analogy  with  the  above  predator-prey  ecological  situation  in 
the  first  model  below.  Qualitative  correspondence  to  the  region-invasion  scenario 
and  simplicity  and  transparency  are  its  major  virtues.  It  is  blatantly  over¬ 
simplified,  a  fault  that  is  partially  rectified  in  Model  II. 

2.1  Mathematical  Model  I:  the  "SCUDWORM"  Problem 

What  follows  is  a  small  modification  of  the  work  of  Ludwig.  Let 
A{t)  =  number  of  attackers  in  !^at  f; 

D{t)  =  number  of  defenders  in  (^at  t. 

Here  t  is  any  real  number,  although  discrete-time  versions  of  the  problem  are 
certainly  possible,  and  possibly  attractive. 

Stipulate  that  these  differential  equations  describe  the  evolution  of  the 
populations: 

A(t)lA)  -A<D(D(t)/l  +  D(()K  (2.1) 

(logistic)  arrival  rate  of  kill  rate  of 

attackers  into  ^  attackers  by  defenders 

pA(t)  -  emim  .  '  (2.2) 

rate  of  rate  of  defender 
defender  force  removal 

force  from  ^ 

allocation 
to  region  ^ 

The  above  setup  cavalierly  attributes  certain  simplistic  behavior  to  both  attacker 
(logistic  increase  in  the  region),  and  the  defender  (attacker  queuing  for 
deterministic  service).  More  realistic  and  inclusive  representations  appear  in  the 
next  model. 
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The  model  of  (2.1)  and  (2.2)  can  be  modified  to  allow  attackers  to  enter  the 
region  when  A{t)  =  0  and  for  numerical  stability  as  follows. 

^  =  [A„  +  M(t)Il  -  A(t)/A]  -  #iD[D(f)/(l + D(f))]pK  (2.2a) 

^  =  M(f)-eD(t)/(A(f)  +  c)  (2.2b) 

where  c  is  a  small  constant. 

The  rightmost  term  of  (2.1)  deterministically  represents  the  expected 
"service"  rate  of  a  D -server /missile  shooter;- see  Filipiak  (1988)  which  relies  on 
the  ideas  of  Rider  (1967)  and  Agnew  (1976).  The  first  right-side  term  of  (2.2) 
states  that  defender  force  allocation  to  the  region  is  enhanced  at  a  rate 
proportional  to  the  current  attacker  population  size;  presumably  refinable  to 
P  =  where  0"  is  a  decision  parameter  and  ^  is  the  rate  of  Blue  sensor  coverage 
of  %  The  second,  or  rightmost,  term  of  (2.2)  postulates  that  defense  force 
depletion  is  proportional  to  that  force,  but  inversely  proportional  to  attacker 
force  size:  if  the  latter  becomes  small,  defenses  become  smaller.  These  rules  are 
arbitrary  but  plausible  (e.g.  0A{t)  could  be  replaced  by  almost  any  increasing 
function  of  A{i),  and  dD{t)/A(t)  by  a  decreasing  function  of  A(t),  without  altering 
qualitative  effects).  It  is  again  acknowledged  that  certain  limited-information  or 
incomplete  battle-damage  assessment  (BDA)  is  not  represented  in  the  present 
abbreviated  model.  For  those  features  see  Model  II,  detailed  in  Section  3. 

Look  for  the  possible  stable  points  of  (2.1),  (2.2).  These  (A  =  A{°°),  D  =  D(oo)) 

‘  » 

satisfy 

0  =  M(i-A/A)-«9kDy^  (2.3) 


O  =  0A-eDIA. 


(2.4) 


Substituting  the  solution  of  (2.4),  i.e.  D  = 


into  (2.3)  we  get 


X4(l->i/I)  =  ff>KD-W  (2.5) 

a  cubic  with  the  possibility  of  one  or  three  real  roots.  Notice  that  if  we  apply  the 
quasi-steady  static  approximation  (QSSA),  cf.  Segel  and  Slemrod  (1989),  letting 
and  6  be  relatively  fast  in  (2.2)  and  putting  the  result  into  (2.1),  the  result  is 


dA{t)  _ 
dt 


A^(t) 


1  + 


iyi‘) 


(2.6) 


which  is  equivalent  to  the  original  model  proposed  by  Ludwig  and  hence  subject 
to  nearly  the  same  analyses,  accessible  in  Murray  (1989).  These  are  based  on 
identifying  A(f)  =  A  =  0  as  one  stable  point  and  factoring  it  out,  then  examining 


i-a/a=|pkD 


(P] 

[ej- 

A 

kJ 

(2.7) 


graphically,  i.e.  plotting  the  left-hand  and  the  right-hand  sides: 
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The  crossings/ intersections  of  left-hand  and  right-hand  sides  of  (2.6)  are  stable  or 
local  equilibrium  points.  The  curve  (a)  can  represent  relatively  rapid  and 
effective  shooter  service,  and  so  provides  for  a  small  equilibrium  attacker 
population  (intersection  of  {£)  and  (a)).  The  curve  (c)  represents  the  opposite: 
slow  shooter  service  or  perhaps  inadequate  sensor  coverage  (P  =  with  ^ 
small)  hence  a  single  equilibrium  point  very  near  the  maximum  regional 
occupancy  chosen  by  Red  (intersection  of  {£)  and  (c)).  Finally,  the  curve  (b) 
represents  the  interesting  bistable  case:  the  basic  parameters  are  consistent  either 
with  a  small  attacker  population  (provided  the  initial  population  is  small),  (left¬ 
most  intersection  of  {£)  and  (b)),  or  a  much  larger  population  if  the  initial 
population  size  is  near  the  rightmost  intersection  of  (i)  with  (b). 

It  is  suggestive  to  "stochasticize"  the  above  system,  replacing  the  system  (2.3) 
and  (2.4),  or  more  accurately  (2.6)  by  a  birth-death  Markov  process  approxi¬ 
mation  that  possesses  a  stationary  distribution.  See  the  appendix  for  details. 

3.  Model  II:  Model  with  Arrivals  and  Departures  of  Attackers  and 

Defenders  and  Simplified  Surveillance:  Delayed  BDA 

Consider  this  more  elaborate  and  realistic  model:  Red  units  (attackers)  enter  ^ 

at  a  specified  rate  which  is  proportional  to  the  number  of  Red  units  that  are  in  the 

#  _ 

region  and  becomes  0  when  the  number  of  Red  units  reaches  a  maximum  A ,  i.e. 
following  logistic  growth  as  in  (2.1);  (modifications  allowed  without  violating 
qualitative  effects).  The  Red  units  are  the  subject  of  a  prescribed  level  of 
surveillance/ reconnaissance,  are  identified  (possibly  incorrectly)  as  alive,  and  are 
eventually  targeted  by  members  of  a  group  of  Defenders  (Shooters).  The 
Defenders  also  arrive  at,  or  focus  assets  on,  the  region  at  a  rate  proportional  to 
the  number  of  attackers  that  have  been  detected.  The  arrival  rate  becomes  0  once 
a  prescribed  number  of  defenders  is  in  the  region. 
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Note  that  the  present  model  explicitly  accounts  for  realistic  effects  that  Model 
I  explicitly  ignores,  namely  imperfect  BDA,  track  loss,  defender  evasion  (leaving 
the  region  to  shake  off  pursuit  -  possibly  inferred  or  imagined  -  and  other  such 
realisms). 

Parameters  for  Model  II 

A  =  Maximum  number  of  attackers  in  the  region 

D  =  Maximum  number  of  defenders  in  the  region  %  This  is  the  number 
of  individual  shooters  that  can  engage  attacker  targets  at  a  time. 

=  Arrival  rate  of  defenders  to  the  region 

6  =  Departure  rate  of  defenders  from  the  region  ^ 

Cs  =  constant 

A  =  Attacker  arrival  rate  into  region  %  Although  this  is  first  treated  as  a 
positive  constant,  it  can  be  made  a  function  of  time,  or  even  of  the 
state  of  Red  forces  in  the  region.  For  example.  Red  may  wish  to 
amass  a  certain  force  size  by  a  particular  time  to  oppose  an 

amphibious  landing.  In  short,  A  again  stands  for  a  simple  control 
variable  ("knob")  that  in  reality  can  be  adjusted  by  Red. 

fi  =  Attacker  attack /"service"  rate  (l//i  =  mean  time  to  track,  shoot,  flight 
time  of  missile). 

V  =  Track-loss  rate  (1/  v  =  mean  of  a  holding  time  of  attacker  in  track).  • 

a  =  Attacker  (e.g.  TEL  =  SCUD  launcher)  shoot  rate  (1/a  =  mean  time 

between  shots  by  a  single  attacker).  This  parameter  is  irrelevant 
under  certain  circumstances;  it  is  relevant  to  SCUD-like  attack. 

Y  =  Rate  attackers  leave  region  (1/7=  mean  holding  time  of  attacker  unit 

in  region).  An  attacker  decision  variable.  . 

^  =  Rate  undetected /tracked  attackers  are  acquired  by  sensors/C^ 

system. 

PK  =  Probability  an  attacked  target  is  killed.  This  parameter  is  currently/ 
presently  taken  to  be  non-range-dependent. 
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Raa  =  Probability  that  an  alive  target  that  is  correctly  classified  as  being 
alive  by  the  sensors /C^  system.  1  -  Raa  =  Rad  =  P(target  classified  or 
perceived  as  dead  |  alive). 

Rda  =  Probability  that  a  dead  target  is  misclassified  as  being  alive  (hence 
potentially  retargeted;  effectively  a  decoy). 

State  Variables  for  Model  II 

The  following  state  variables  are  needed  to  describe  the  above  dynamical 
system; 

Au(t)  =  Number  of  undetected  live  (hence  potentially  active  and  threatening) 
Red  attackers  present  in  region  <^at  time  t. 

Ad(t)  =  Number  of  detected  live  Red  attackers  present  at  time  t.  These  are  on 
the  Blue  shooter's  target  list,  and  will  be  engaged  unless  lost  by  the 
sensor  system  (they  may  go  into  hiding,  or  even  leave  the  region 
covered  by  the  surveillance,  e.g.  JSTARS). 

Da(t)  =  Number  of  detected  and  perceived  to  be  alive,  hence  potentially  Blue- 
targeted,  but  actually  dead  Red  attackers  at  f.  These  are  present 
because  Blue  battle  damage  assessment  (BDA)  is  realistically . 
imperfect. 

Dy(f)  =  Number  of  dead  attackers  in  region  that  are  not  yet  classified. 
Classification  to  be  done  by  the  surveillance/reconnaissance  system. 

S(t)  =  Number  of  defenders  (shooters)  in  the  region  ^ 


For  the  deterministic  modeling  of  saturable'  service  by  the  defense 
shooters  (presumably  naissile  launchers  in  the  present  context)  we  make  use  of 
the  following  approximation,  cf.  Filipiak  (1988),  Agnew  (1976),  and  also  Rider 
(1967),  for  the  saturable  rate  of  processing  by  the  shooters, 

"s'*'-  1+[A,(()  +  D.(t)] 
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Transition  Equations  for  Model  II 


2  + 


A 


+vA^(f)  -  Au(t}(a  +  r)  -  ^RaaM^) 

=  ^oaM^)  +  «^«(0  -  (r  +  v)Arf(f) 


^Ml- 

dt 


^a{t) 


^  =  PK(f)  +  D„(f))[l-§i 


-0 


5(0 


^rf(0  "*"  ^fl(0  ^ 


(3.1a) 


(3.1b) 


(3.1c) 


,  '«s(<)Pk  -|D|.(t)+  .■  /  .Hs(t)  (3-ld) 


(3.1e) 


Solutions 

It  is  possible  to  obtain  some  initial  information  in  the  form  of  steady-state  or 
long-run  solutions.  Set  the  derivatives  of  the  equations  (3.1a)  -  (3.1e)  equal  to  0  to 
find: 


0  =  A[A„  +  A^] 
Ad 


2  ^  A^ 


+— ^Hs[l-px]  +  vA^-A„(a  +  r  +  ^J 


Arf  +l^fl 


(3.2a) 


0  -  -{Y  +  v)A^  -  --— . . Hg 

Ad 

0  =  «R*D„-— 

■^d  ■*■ 


(3.2b) 

(3.2c) 
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0  = 


0  =  P(A^  +  D,) 


1-1' 

D 


-e 


+Dn  +  Cc 


M  '  ^a  '  ^s 

Solving  (3.2d)  for  Du  and  substituting  the  expression  into  (3.2c)  results  in 


D„=i 


Ad  --  O'’ 


“  ^Ud+D, 


-^sPk  + 


Ad  +Da 


-He 


and 


0  =  R, 


da" 


Aa  „  D, 


D„ 


HsPK+-^^Hc:\—^^He 


Aa+Du~~^^-  A^+Du  "J  + 
Ofl  (l  -Rda)  =  ^da^dPK 


1-R. 


■da 


Sdlving  (3.2e)  for  S  results  in 


S  = 


P{A^  +  Da)  _  ^(1  +  c)A^D  ((1  +  c)A^  +  Cg ) 

^(A^Dg)  Q  P[l  +  c)Ai}{(l  +  c)A^-i-Cg^  +  OD 

D  Aa+D^+Cs 


Equation  (3.2b)  results  in 


^u=- 


1 


Ad 


(y  +  v)A^  +  — ; — —  piS 

U  )  d  i  +  (i  +  c)A^^  ^ 


(3.2d) 

(3.2e) 


^aa  "*■  ^  1 

Substitute  Au  into  equation  (3.2a)  and  (numerically)  solve  for  Ad- 
Numerical  Example 

For  the  parameters  Rad  =  Rda  =  0.3,  =  10,  A  =  10,  D  =  10,  ^  =  500,  a  =5, 

PK  =  0.7,  V  =  0.5,  Y=  0.1,  6  =  10,  Cs  =  0.1,  p  =  20,  the  steady  state  equations  (3.2a)  - 
(3.2e)  have  4  solutions 

{Au/  Ad,  Da,  Du,  S)  =  (0, 0, 0, 0, 0) 

{Au,Ad,Da,Du,S)  =  (0.04,1.07,0.32,0.03,1.72) 
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{AurAd,Da,Du,S)  =  (0.4,13.4,4.0,0.28,9.7) 

{AurAci,Da,Du.S)  =  (0.48,33.7,10.1,0.30,9.95). 

The  solutions  (0,  0,  0,  0,  0)  and  (0.4,  13.4,  4.0,  0.28,  9.7)  are  unstable 
equilibrium  points.  The  solutions  (0.04, 1.07,  0.32,  0.03, 1.72)  and  (0.48, 33.7, 10.1, 
0.30, 9.95)  are  stable  equilibrium  points.  That  is,  for  the  present  example  there  are 
two  distinct  steady  states:  irdiabitation  of  one  or  the  other  depends  on  where  the 
system  started  from,  i.e.  the  initial  conditions. 

Thus,  the  differential  equations  can  have  different  limiting  solutions 
depending  on  the  initial  conditions,  a  classical  multistable  situation.  In  particular, 
for  the  present  numbers: 

if  A„(0)  =  5,  AdiO)  =  0,D„(0)  =  0,  S(0)  =  0,  D„(0)  =  0, 

then  AuH  =  0.05,  AdH  =  1.07,  DaH  =  0.32,  SH  =  1.73;  D„(oo)  =  0.031. 

In  this  case  the  defenders  were  able  to  oppose  the  attackers  before  buildup 
occurred,  and  to  force  the  number  present  to  nearly  zero.  This  might  well  be  the 
fruit  of  superior  intelligence  and  the  available  assets  to  capitalize  on  it.  However, 
if  A„(0)  =  30,  AdiO)  =  0,  DaiO)  =  0,  S(0)  =  0,  D„(0)  =  0, 
then  AuH  =  0.48,  AdH  =  33.7,  DaH  =  10.1,  SH  =  9.95,  DuH  =  0.30. 

Then  in  this  case  the  Blue  force  was  presumably  oblivious  to  attackers  until  their 
number  became  substantial,  (possibly  a  failure  of  intelligence  and/ or  overall 
surveillance)  after  which  defenders  swung  into  action  and  detected  and  targeted 
those  defenders  in  the  region.  The  number  of  shooters  appears  inadequate 
because  the  number  of  attackers  queued  for  shooter  service  eventually  is  higher 
than  the  initial  number  in  the  region.  If  assets  are  available  this  deficiency  might 
well  be  corrected  in  real  time.  But  the  present  model  allows  anticipation  of  such  a 
possible  state  of  affairs,  and  might  allow  forestalling  it.  This  is  "real  dynamics" 
that  is  unmodeled,  but  for  which  models  are  useful. 


12 


REFERENCES 


Agnew,  C.E.  (1976).  Dynamic  modelling  and  control  of  congestion  prone 
systems.  Operations  Research,  24,  pp.  400-419. 

Dockery,  J.T.  and  Woodcock,  A.E.R.  (1993).  The  Military  Landscape:  Mathematical 
Models  of  Combat.  Woodland  Publishing  Ltd.,  Cambridge,  England. 

Filipiak,  J.  (1988).  Modeling  and  Control  of  Dynamic  Flows  in  Communication 
Networks.  Springer-Verlag,  Berlin. 

Gaver,  D.P.  and  Jacobs,  P.A.,  A  model  for  analyzing  Blue  force  response  to  region 
invasion  by  multi-type  Red  forces.  Naval  Postgraduate  School  Technical 
Report,  forthcoming. 

Gaver,  D.P.  and  Jacobs,  P.A.  (1998).  Enemy  air  defense  (EAD)  suppression  as  an 
information  duel.  Naval  Postgraduate  School  Technical  Report  NPS-OR-98- 
005. 

Ilachinski,  A.  (1996).  Land  Warfare  and  Complexity,  Part  II:  An  Assessment  of  the 
Applicability  of  Nonlinear  Dynamic  and  Complex  Systems  Theory  to  the  Study  of 
■Land  Warfare..  Center  for  Naval  Analyses,  Alexandria,  VA. 

Ludwig,  D.,  Aronson,  D.G.,  and  Weinberger,  H.F.  (1979).  Spatial  patterning  of 
the  spruce  budworm.  J.  Math  Biology,  8,  pp.  217-258. 

Ludwig,  D.,  Jones,  D.D.,  and  Rolling,  C.S.  (1978).  Qualitative  analysis  of  insect 
outbreak  systems:  the  spruce  budworm  and  forest,  f.  Anim.  Ecology,  47,  pp. 
315-332. 

Murray,  J.D.  (1989).  Mathematical  Biology.  Springer-Verlag.^ 

Rider,  K.L.  (1967).  A  simple  approximation  to  the  average  queue  size  in  the  time- 
dependent  queue.  Journal  of  the  ACM,  23,  pp.  361-367. 

Segel,  L.A.  and  Slemrod,  M.  (1989).  The  quasi-steady  state  assumption:  case 
study  in  perturbation.  SIAM  Rev.,  31,  pp.  446-447. 


13 


APPENDIX 

A  Birth-Death  Markov  "SCUDWORM"  Model 
Suppose  that  A(f)  now  represents  a  Markov  chain  in  continuous  time  with 
state  space  {0,1,2, ...,  A}.  Here  A  is  an  integer.  Furthermore,  let  defender 
adaptation  dynamics,  as  represented  by  P  and  0  be  fast  as  compared  to 
parameters  A  and  jn,  and  invoke  QSSA.  Then  we  replace  the  deterministic 
differential  equation  (2.6)  by  a  birth-death  process  with  generator 

P{A{t  +  dt)  =  A{t)  +  l|A(f)}  =  AA(f)(l  -  A{t)/A)dt  +  o{dt)  (A.l) 

and 

Di^]A^{t) 

P{A{t  +  dt)  =  A{t)  - 1| A(f )}  =  upk  - ■  dt  +  o{dt) ;  (A.2)- 

there  is  a  natural  boundary  at  A(f)  =  0  so  no  transition  to  negative  values  is 
possible.  Additionally, 

r  _  D(|u^(f) 

P{A(t  +  dt)  =  A(t'iA(t)}  =  l-  \A{t){l-A(t)/A)  +  tv>K  ■  dt  +  o{dt).  (A.3) 

All  other  transitions  have  negligible  probability  o(dt)  as  0.  Because  of  the 
non-linearity  of  the  transition  functions  it  is  not  true  that  £[A(t)]  satisfies  the 
differential  equation  (2.6),  although  there  will  be  qualitative  resemblance 
between  E[A{t)]  and  the  differential  equation  solution  A{t)  as  functions  of  f.  It 
should  also  be  possible  to  show  that  A(f)/  A  tends  weakly  to  a{t)  =  A{t)/  A  as  A 
becomes  large,  provided  D  =  d  A;  furthermore  it  should  be  the  case  that,  for 
large  A,  A(f)  can  be  well  approximated  as  the  sum  of  a  deterministic  trajectory 
and  diffusion  process.  This  option  is  not  investigated  here;  something  simpler  is. 
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The  stochastic  model  (A.l),  (■A.2),  (A.3)  is  clearly  richer  than  its  deterministic 
counterpart  or  approximation,  (2.1)  -  (2.6).  Furthermore,  the  present  simplified 
stochastic  version  can  be  studied  reasoriably  explicitly  and  analytically,  i.e. 
without  Monte  Carlo  simulation  or  intricate  numerical  computation.  Various 
questions  will  now  be  addressed.  For  simplicity  put 

P|A(t  +  dt)  =  a  +  l|A(t)  =  fl}  =  Xgdt  +  o{dt)  (A.4) 

and 

P|A(f  -\-dt)  =  a- 1|  A(f)  =  fl}  =  fi^dt  +  o{dt)  (A.5) 

where  Xa  and  ya  are  gotten  from  (A.l)  and  (A.2).  Look  at  some  options. 

(1)  Stationary  solution 

As  (A.l)  is  presented,  it  describes  an  absorbing  chain  with  certain  absorption 
at  A(f)  =  0.  Modify  this  by  introducing  an  attacker  import  rate  Aq  (Aq  >  0):  with 
probability  Xodt  a  new  attacker  appears  if  the  region  becomes  empty.  Now  using 
well-known  formulas  one  can  find 

lim  P{  A{t)  =  A(0)  j  =  TtQ  =  7ra,  0<a<  A,  (A.6) 

f— ><» 

where  tcq  normalizes  iia  to  sum  to  unity.  This  probability  mass  fimction  Tia  can  be 
studied  analytically,  or,  more  easily,  evaluated  and  studied  numerically. 
Exploration  shows  that  it  can  be  bimodal,  with  modes  corresponding  roughly  to 
the  two  stable  points  that  may  occur  for  the  deterministic  model.  An  example  is 
provided  in  Figure  Al.  The  operational  interpretation  is  that,  over  a  long  period 
of  .time  attacker  numbers  may  typically  reside  near  either  a  low  value,  or  a 
higher  -  obviously  more  dangerous  -  level.  The  chance  of  being  at  that  level  can 
be  controlled  by  enhancing  defense  shooter  assets.  In  a  rough  sense  the  A(f) 
process  tends  to  be  near  one  or  the  other  of  the  two  stable  points,  remaining  at 
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one  place  for  a  nearly  exponential  time  period,  then  jumping  to  (near)  the  other, 
and  so  on  back  and  forth  until  the  parameters  change.  The  defender's  job  is  to 
cost-effectively  lengthen  the  time  attacker  state  resides  at  a  low  level,  while  the 
attacker  presumably  wants  her  state  to  remain  high,  but  also  to  carry  out 
aggressive  acts,  such  as  launching  SCUDs. 

(2)  First-Passage  or  Exit  Times 

Suppose  that  initially,  i.e.  at  f  =  f'  (where  t'  can  be  taken  to  be  0  because  of 
stationary  transitions)  A{t')  =  i,  any  particular  state  value.  Then  it  is  of  interest  to 
study 

tij  =  inf{f :  A{t  + 1')  =  j\A{t')  =  /}  (A.7) 

the  first-passage  time  (fpt)  from  i  to  j  ^  i;  this  is  the  random  time  required  to  pass, 
for  the  first  time,  from  state  A(f')  =  i  to  state  j.  For  models  of  the  present  type  it  is 
well-understood  how  to  calculate  all  interesting  probabilistic  features  of  Ty.  For 
example, 

(2.1)  Expected  fpt.  Let  mij  =  £[Ty].  Then  it  is  easily  shown  that  my  satisfies  a 
backward  equation: 


with  boundary  condition  m/f  =  0,  mgj  =  0  if  ^  <  0. 

Example:  Suppose  A{t')  =  0.  Find  £[to2]  =  wo2- 
From  (A.8), 

mo2  =  q— - — [^12^1  +  m_i^2A^-i  + 1] 
Aq+Pq'- 

=  ^[mi2Ai+l] 

by  virtue  of  boundary  conditions.  Also 


(A.8) 


(A.9.1) 
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mi2  =  - - [^22^2  +  mQ2fli  + 1] 

Ai  +  jUl 


(A.9.2) 


1 

Ai  +  /zi 


[mo2iUi+l]. 


The  two  equations  (A.9.1)  and  (A.9.2)  can  now  be  solved  simultaneously.  This 
pattern  can  be  followed  more  generally. 

Such  stochastic  models  are  usually  not  available  in  a  form  even  as  simple  as 
the  above.  Monte  Carlo  simulation  will  ordinarily  be  required  in  order  to  elicit 
lessons  and  discover  sensitivities.  But  a  quite  abbreviated  model  such  as  the 
present  "SCUDWORM"  example  is  useful  in  that  it  may  quickly  indicate 
particular  sensitivities  and  tradeoffs,  in  this  case  to  increased  surveillance 
defense  surveillance  or  increased  attacker  activity  (perhaps  putting  =  (|  +  oc)^, 
where  ^  is  surveillance  rate  and  a  is  attacker  rate  of  fire). 
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NUMBER  OF  ATTACKERS 


LIMITING  PROBABILITY 
0.04  0.08  0.12 
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